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^ i Abstract 

| In the causal theory of relativistic dissipative fluid dynamics, there are conditions on the equation 

o : 

of state and other thermodynamic properties such as the second-order coefficients of a fluid that 
need to be satisfied to guarantee that the fluid perturbations propagate causally and obey hyper- 

O 

bolic equations. The second-order coefficients in the causal theory, which are the relaxation times 



for the dissipative degrees of freedom and coupling constants between different forms of dissipation 



(relaxation lengths), are presented for partonic and hadronic systems. These coefficients involves 
relativistic thermodynamic integrals. The integrals are presented for general case and also for dif- 

•i-H . 

ferent regimes in the temperature-chemical potential plane. It is shown that for a given equation 
of state these second-order coefficients are not additional parameters but they are determined by 
the equation of state. We also present the prescription on the calculation of the freeze-out particle 
spectra from the dynamics of relativistic non-ideal fluids. 

PACS numbers: 05.70.Ln, 24.10.Nz, 25.75.-k, 47.75. +f 
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I. INTRODUCTION 



Transport properties of relativistic hot and dense matter are important in many physi- 
cal situations such as in nuclear physics, astrophysics and cosmology. The hot and dense 
relativistic matter is produced in high energy relativistic nuclear reactions such as those in 
heavy-ion experiments while the hot matter is believed to have existed in the early universe 
and the dense matter might be found in the quark stars and neutron stars. The transport 
properties of a system out of equilibrium are governed by transport coefficients which re- 
late flows to thermodynamic forces. These coefficients characterize the magnitude of the 
response of the system (flow) to a certain disturbance (thermodynamic force). 

The space-time evolution of the produced matter in high energy nuclear collisions can 
be studied using relativistic fluid dynamics. In order to take into account non-equilibrium 
(dissipative) effects we need to use relativistic non-ideal fluid dynamics. We consider the 
causal theory of relativistic fluids as presented in paper I [1]. 

In the high energy heavy ion collisions the rapid evolution of the fluid is governed by 
dissipative effects via transport coefficients such as the shear and bulk viscosities, the ther- 
mal conductivity, the diffusion coefficients, etc. These transport coefficients are generally 
calculated via the use of kinetic theory and thereby imply the knowledge of a collision term. 
The collision term must be a representation of the hot and dense matter at hand. 

In addition to the standard transport coefficients the causal theory contains other thermo- 
dynamic functions, namely the second-order coefficients. These coefficients, in combination 
with the standard transport coefficients, are related to the relaxation times and relaxation 
length for various dissipative processes. The relaxation length results from the coupling 
between the heat flux and viscous processes. These second-order coefficients are given by 
complicated thermodynamic integrals. In this work we will present the second-order trans- 
port coefficients in detail for future applications in the description of the relativistic dynamics 
of heat-conducting and viscous matter produced in high energy nuclear collisions. 

In addition to the knowledge of the transport coefficients and second-order thermody- 
namic function we need to know how to determine the form of the equation of state. Knowl- 
edge of the equation of state is needed for analyzing many important physical situations such 
as supernova explosions and neutron star formation, high energy nuclear collisions, and the 
study of quark-hadron phase transition in the early universe. Recent developments in heavy 
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ion collision experiments can reveal the form of the equation of state of the nuclear matter. 
The physics involved in such collisions is highly complex and it is highly nontrivial matter 
to extract equation of state information from such a complex dynamical situation. 

In this paper we would like to investigate the influence of the equation of state on the 
transport properties of the relativistic hot and dense matter. The equation of state in the 
hadronic phase is derived from the Walecka model of nuclear matter and in the QGP phase 
we use the MIT Bag model. In the mixed phase we employ the Gibbs construction for the 
phase equilibria. 

The rest of the paper is organized as follows: In Section II we present the moments of 
relativistic Boltzmann equation in terms of the relativistic thermodynamic integrals. In Sec- 
tion III we present the 14-Field Theory from the microscopic considerations. We consider 
small deviations from equilibrium and employ the Grad's 14 moment method. In Section 
IV we present the second order entropy 4-current in terms of the relaxation and coupling 
coefficients which are collectively referred to as the second order coefficients. In Section V 
we present the EoS considered in this work. In Section VI we present the freeze-out pre- 
scription which takes into account the dissipative corrections. In Appendix A we present the 
relativistic thermodynamic integrals. In Appendix B we present the various limiting cases 
of the thermodynamic integrals. In Appendix C we present the ultra-relativistic thermody- 
namic integrals. In Appendix D we present the relativistic thermodynamic integrals at zero 
temperature. 

Our units are h = c = k B = I- The metric tensor is g^ u = diag(+, — , — , — ). The scalar 
product of two 4- vectors a M , 6 M is denoted by a^g^ u b u = a M 6 M , and the scalar product 
of two 3-vectors a and b by a ■ b. The notations A {ap) = (A af3 + A /3a )/2 and A [a ® = 
(A a @ — A f5a )/2 denote symmetrization and antisymmetrization, respectively. The notation 
A ( af 3) = [ A (a A /3) A «/3 A^]A^ denotes the trace-free part of A^ . The 4-derivative is 
denoted by d a = d/dx a . An overdot denotes A = u x d\A. 

II. THE RELATIVISTIC TRANSPORT EQUATION AND THE MOMENT 
EQUATIONS 

The transport equation gives the rate of change of the distribution function in time and 
in space due to the particle interactions. The distribution function can in principle be solved 
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from the transport equation. Here we limit ourselves to cases valid for dilute systems where 
binary collisions dominate. In trying to study high energy nuclear reactions we need to know 
the properties of many-particle system. While the properties of such a system depends upon 
the interactions of the constituent particles and external constraints in kinetic theory in the 
macroscopic level such a system is described by the net conserved densities, the energy 
density, hydrodynamic velocity and dissipative quantities. Thus the resulting fluid dynamic 
equations can be considered as an effective kinetic theory. The results should be compared 
with other theories which goes beyond dilute systems to check the deviations. 

The relativistic Boltzmann transport equation for the invariant on-shell phase space den- 
sity f(x,p) is 

P^dJ{x, p) = jf Ult Df(x, p) + p^ Vm f{x, p) = C[f(x, p)\ , (2.1) 

where C[f] stands for the collision term. Here D = w M <9 M is the covariant time derivative 
and — A^ct is the covariant gradient operator. Furthermore, A Miy = g^ v — u^u" is 
the projection tensor and w M = (7, 71;) is the 4-velocity where 7 = y/l — v 2 and the three- 
velocity of a particle is v = p/p°. The four velocity is normalized such that u^u^ = 1. The 
four momentum is p^ = (p°, p) where p° = E = \/p 2 + m 2 is the relativistic energy of the 
particle and m is the mass of the particle. Eq. (2.1) describes the time evolution of the 
single particle distribution function f(x,p). The local equilibrium distribution function has 
the form 

n(^P)=Ao e ^ 1 _ a _ a (2.2) 

where A = g/(2-rc) 3 and a corresponds to the statistics of Boltzmann (a = 0), Bose (a = +1), 
and Fermi (a = —1) distributions. Also the degeneracy g = (2 J + 1) where J is the spin of 
the particle, (3^ = flu^, a = fl/i with (5 = 1/T the inverse temperature, and \i is the chemical 
potential. 
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Let us define the first five moments of the distribution function f(x,p) by 

N»(x) = J f(x,p)jfdw , (2.3) 

T^(x) = J f(x,p)p^p l/ dw , (2.4) 

F Xpu (x) = J f(x,p)p x p p p u dw , (2.5) 

R a ^ v {x) = J f{x ) p)p a p^p^p u dw , (2.6) 

Q a ^ up {x) = J f{x,p)p a p^p^p u p p dw , (2.7) 



where 

dw = (2.8) 

For the equilibrium distribution function f cq (x,p), Eq. (2.2), the five moments can be written 

as 

= Ziou" , (2.9) 

T% = l 20 u>*u»-l 21 A^ , (2.10) 

F Xp » = X m u x u p u v - 3J 31 u v) , (2.11) 

= Z m u a u p u»u v - 6J 4 i A (al3 u p u u) + 3J 42 A (a/3 A pv) (2.12) 

Qcfavp = x 50 u a u^u u u p - IOJ51 A^ vrvrvfl + 15J 52 A^u^ (2.13) 

where the T nk are the equilibrium thermodynamic functions and are presented in Appendix 
A. Various contractions of the above equations produce useful relations such as 

=I 10 = n, u^u v T% = X 20 = e, T ccl p = J 20 - 3J 21 = e - 3p , (2.14) 
u x F eq Xv = J 30 - 3J 31 = m 2 n , u a u p R C(l f u = J 40 - 3J 41 = m 2 e . (2.15) 

Corresponding to the five moments are the five auxiliary moments of A(x,p)f(x,p), which 
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arises due to variations in the distribution function, 

A>(x) = J A(x,p)f(x,p)p**dw , (2.16) 

f^(x) = J A(x,p)f(x,p)p^dw , (2.17) 

F A ^(x) ee y A(x,p)/(i,p)pV/d«; , (2.18) 

r*Pi»>( x } = J A(x,p)f(x,p)p a p l3 p^p u dw . (2.19) 

q°0i»>p( x ) = J A(x,p)f(x,p)p a pPp»p u p p dw , (2.20) 

where A(x,p) = 1 + aAp 1 f(x,p). For the equilibrium distribution function, Eq. (2.2), and 
A eq (x,p)f eq (x,p) the auxiliary moments can be written as 

AT q = Jiou" , (2-21) 

= t 7 20 u'V - J21 A"" , (2.22) 

= J 30 u x uV - 3 J 31 A< A " u v) , (2.23) 

R«P^ = J^ u a u p u^u v - 6J41 A (a/3 + 3 J" 42 A (a/3 A""> (2.24) 

ga^p = j bGU * u P u » u " u P _ ioj- 51 A( a/3 u^u") + 15 J52 A( a/3 . (2.25) 

where the relativistic thermodynamic integrals J n k are also presented in Appendix A. Note 
that by Eq. (A12) the J n k can be related to the X n k, for example 

rr _ ZlO rr (^20 + ^2l) ~ _(J 3 0+2J 3 l) _J 3 1 

J21 -y, J31- ~ p , J4i- ^ , J42 -y (2.26) 

Also by Eq. (A13) the J n k can be obtained as differentiation of the X n k with respect to a 
and (3. 

As in the case of the X nk similar contractions of the J n k leads to similar relations, Eqs. 
(2.14)-(2.15), with the T nk now replaced by the J n k- For later reference we define the entropy 
4-current in kinetic theory as 

S»(x) = - J dwjf {f(x,p) In (A Q l f{x,p)) - a - l A Q A{x,p) In A(x,p)} , (2.27) 

and we introduce a quantity 

h, = £±£> = » = , (2.28) 

™i Tlj J 2 1 

where w = 5 + p, as the enthalpy per net conserved ith charge. Since in our case we are 
considering one such charge, namely the net baryon number, h denotes the enthalpy per net 
baryon. Similarly one defines s/n as the entropy per net baryon. 
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III. THE 14-FIELD THEORY OF NON-EQUILIBRIUM RELATIVISTIC FLUID 
DYNAMICS 

A. Small deviations from thermal equilibrium 

For a gas that departs slightly from the local thermal equilibrium, we may write the 
distribution function as 

f(x,p) = r(x,p) {1 + A cq (x,p)(f)(x,p)} , (3.1) 

where 4>(x,p) is the deviation function to be discussed in the next subsection. Substitution 
of Eq. (3.1) into Eqs. (2.3), (2.4), and (2.27) yields 

N»(x) = N^(x) + 6N»(x) , (3.2) 
T^(x) = T%{x) + ST^ix) , (3.3) 
S"(x) = S^(x) + 5S»(x) , (3.4) 

where Ng q , Tg£, and Sg q have the ideal fluid (equilibrium) form, and 

5N"(x) = J / eq (i,p)A cq (i,j))^i,ri/dw , (3.5) 

8T^{x) = J f cq (x,p)A cq (x,p)<j ) (x,p)p> 1 p u dw . (3.6) 

For the entropy, expanding the term in the curly brackets in Eq. (2.27) up to terms of second 
order in deviations 4>(x,p) (cf. Section IV) leads to 

5S"(x) = S"(x) - S^x) 

= - J [a(x) - p v {x)p»] r(x,p)A cq (x,p) ( l ) (x,p)p^dw 

-\ J f cq (x,p)A cq (x,p) (<f>(x,p)) 2 p^dw + • • • (3.7) 
which corresponds to the phenomenological expression 

5S"(x) = -a{x)5N^{x) + (3 u (x)5T^(x) + Q»(x) , (3.8) 

with 

Q^{%) = ~\j r(x,p)A^(x,p)) 2 I fdw + ■■■. (3.9) 



now explicitly denned in kinetic theory. The five parameters a and /3 M which describe the 
nearby equilibrium state can be specified by matching the equilibrium state to the actual 
state. This is done by the prescription that the net charge density and energy density are 
completely determined by the equilibrium distribution function. That is, we impose the 
following conditions: 

n = J dwp IJ 'u^f(x,p) = n cq = J dwp^u^f^ix, p) , (3.10) 

S EE J dw(p»U,) 2 f(x, P ) — E C q EE J dw ( P » U ,) 2 T (x , P ) . (3.11) 

Then p(e, n) = p(e eq , n cq ) is the actual equation of state. The above matching conditions 
also implies 

UpSNi 1 = u^uJT^ = . (3.12) 
It then follows from Eqs. (3.8) and (3.9) that 

vss** = V?" < , ( 3 - 13 ) 

confirming that equilibrium maximizes the entropy density under the above constraints, 
Eq. (3.12), on particle and energy densities , and that Sg q and Tj^ give the entropy and 
thermodynamical pressure correctly to first order in deviations. To fix the hydrodynamic 
velocity we may choose either the Eckart's or Landau-Lifshitz' definition of the 4-velocity. 
The former case implies 

A" U N U = J dw A^^r^p) = , (3.14) 

and the latter case implies 

AT w u ff = J dw p vPa u a f eq (x,p) = . (3.15) 

In finding the equations for non-ideal fluid dynamics from microscopic/kinetic the- 
ory there are basically two approaches which lead to linearized transport equations: the 
Chapman-Enskog approximation [2] and Grad's 14-moment approximation [3]. In the 
Chapman-Enskog method one solves a linearized Boltzmann equation for f(x,p) which is 
linearized by assuming that the function <f>(x,p) that measures deviation from local equi- 
librium is small. Terms which are non-linear in <f>(x,p) and also the relative variation over 
a mean free path are neglected. One then uses the conservation laws d^N^ = d^T^ = 



to express the time derivatives D of the thermodynamical variables and the four-velocity 
(hence D/ eq (x,p)) in terms of spatial gradients Vm of a ( x ) an d P^(x), correct to first order. 
Substitution into the kinetic expressions (cf. Eq. (2.1)) connects these gradients to heat flow, 
diffusion and viscosity. The second approximation which is based on the moment method 
is more general and has a wide scope of applicability. It will be discussed in the following 
section. 

B. Grad's 14 Moment Approximation 

In the phenomenological theories of relativistic non-ideal fluid dynamics by Eckart [4] and 
by Landau and Lifshitz [5] (see [1] for details] , instantaneous propagation of heat and viscous 
signals (acausality problem) remained a puzzle for many years. However, in 1966, I. Miiller 
[6] traced the origin of the difficulty to the neglect of terms of second order in heat flux 
and viscosity in the conventional theory's expression for the entropy. Restoring these terms, 
Miiller was led to a modified system of phenomenological equations which was consistent 
with the linearized form of Grad's kinetic equations. Midler's theory was rediscovered and 
extended to relativistic fluids by Israel and Stewart [7] in the 1970's. 

The progress made in phenomenology in getting rid of the acausality problem can also 
be made from kinetic theory. The analogous paradox in non-relativistic kinetic theory was 
resolved by Grad [3], who showed in 1949 how transient effects could be treated by em- 
ploying the method of moments instead of the Chapman-Enskog normal solution. Suitable 
truncation of the moment equations gave a closed system of differential equations which 
turned out to be hyperbolic, with propagation speeds of the order of the speed of sound. 
The relativistic version of Grad's method was developed by Israel and Stewart [7]. 

In phenomenology we seek a truncated hydrodynamical linearized description of small 
departures from equilibrium in which only the 14 variables iV M (x) and T^ v (x) appear. The 
microscopic counterpart of this is a truncated description in which the function 



y(x,p) = In [A 1 f(x,p)/A(x,p)] 



(3.16) 



differs from any nearby local equilibrium value 



(3.17) 
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by a function of momenta specifiable by 14 dynamic variables. The truncated description 
is accomplished by postulating the relativistic Grad's 14-moment approximation [3, 7], or 
variational method [8], that y — y eq can be approximated by a quadratic function in momenta 

4>{x,p) = y(x,p) - y ecL (x,p) = e{x) - e^p* + e^x^p" + ... , (3.18) 

or 

y(x,p) = (a + e(x)) - + e^(x))jf + e^{x)p»p v + ... , (3.19) 

where e(x), e p (x), and t pv [x) are small. Without loss of generality e pu (x) may be assumed 
traceless, since its trace can be absorbed in redefinition of e(x). The non-equilibrium distri- 
bution function is given by Eq. (3.1) and it depends on the 14 variables a + e , f3 p + e p , and 
e pu . While five of these determine the equilibrium state the other nine variables are related 
to the dissipative fluxes. Inserting the expression for f(x,p), Eq. (3.1), into the kinetic 
expressions for and T pv ', Eqs. (2.3) and (2.4), we then have 

N" = N^ + e J dwf q (x,p)A cq (x,p)p" 

-e v J dwf cq (x, p) A cq (x, p)// + t vX J dwf eq (x, p) A cq (x, p)p p p v p x (3.20) 
= Tg + e J dwf cq (x,p)A cq (x,p)p> l p u - e x J dwf cq (x,p)A cq (x,p)p»p l/ p x 

+e Xp J dwf cq (x,p)A cq (x,p)p^p l/ p x p f> (3.21) 
p^v = pov» + e J dwf eq {x,p)A cq {x,p)p a p p p u -e x J dwf cq (x,p)A cq (x,p)p a p»p u p x 

+e Xp J dwf cq {x,p)A cq {x,p)p t7 p p p v p x p p (3.22) 

Then by Eqs. (2.16)-(2.19) one can write the non-equilibrium part of the number four current, 
the energy-momentum tensor and the fluxes as 

SN» = eN^ - e x f£ + e Xu F x ^ , (3.23) 
5T^ = - e x F x ^ + e Xp R x ^ , (3.24) 
SF *v = e p*v> _ ex R^ x + e Xp Ql^ xp . (3.25) 

From the definitions of the dissipative fluxes, namely the particle drift flux = A P 5N", 
the energy flux = A^u p 5T" p , the heat flux g M = — hV^, the bulk viscous pressure 
II = — I A pu 5T^ U and the shear stress tensor n pu = ST^ and the matching conditions, 
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Eq. (3.12), one then obtains the 14 variables a+e , /^+e M , and e^ v in terms of the macroscopic 
fields n, e, n*, IT, g", and tt"" [7] 

e M z, = -4 2 (3« M ^ - A^)I1 - B x u {il q v) + C tt^ , (3.26) 
e M = - , (3.27) 

e = An . (3.28) 

The coefficients At, Bi and C are thermodynamic functions given by 



A = SA^o' {4(JioJ4i - J20J31)} , i3 = Bi^i , (3.30) 
A = 3^2^ Pao + J41 J20 - J30J31} , (3.31) 



with 



A = % , (3-32) 



u = , , d\nl 31 \ +5 
d^l 10 J s/n 



= "3 * ^r-== * '- + 5 , (3.33) 

D n k = Jn+\,kJn-\,k ~ (Jnk) 2 ■ (3.34) 

Once the deviation function is determined one can then derive the equations of motion 
for n, e, IT, g M , ir^". In kinetic theory one uses the moments of Boltzmann transport 
equation, Eq. (2.1). For a general tensorial function of momenta we have 

J dw <$>{p)p»d^f{x,p) = J dw $(p)C[f] . (3.35) 

For Grad's 14-moments approximation we only need the equations for the first three mo- 
ments of the distribution function f(x,p). For $(p) = 1, p u , and p u p x we have respectively 

j dw p»dj(x,p) = = J dw C[f] = , (3.36) 

J dw p^dj(x,p) = d,T^ = Jdw p v C\j\ = , (3.37) 
J dw p"p u p%f(x,p) = d^ uX = J dw p u p x C[f] = P uX . (3.38) 
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The first two moment equations give the 5 conservation laws, the particle number and 
energy-momentum conservation. The remaining additional 9 equations are obtained from 
the third moment equation which represents the balance of fluxes. F^ vX is a completely 
symmetric tensor of fluxes and P" u is its production density which is the rate of production 
per unit 4- volume of due to collisions. The balance equations are 15 in number for only 
14 fields. The trace of F" vX is mass squared times Eq. (3.36) and P^ is traceless. That is, 

P v v = and = m 2 N» , (3.39) 

so that the trace of the tensor equation, Eq. (3.38), reduces to the conservation law of particle 
numbers, Eq. (3.36). Among the 10 equations, Eqs. (3.38), there are only 9 independent ones. 
The resulting relaxational transport equations [7] are presented in detail for applications in 
[1]. Thus the set of equations, Eqs. (3.36)-(3.38), is a set of 14 independent equations 
for 14 fields. In this 14-Field Theory (a field theory of the 14 fields) for non-equilibrium 
relativistic fluid dynamics the dynamical equations, Eqs. (3.36)-(3.38), with representations 
(tensor decomposition) in the Eckart or particle frame 

N» = nu", (3.40) 
= eu»u v - (p(e, n) + II) A"" +2q^u v) + (3.41) 

pw = \c n A 2 (3u'V - A"") II + 2C q B iq K u) + \c^C^ v \ (3.42) 
= (ls + [A o J 30 -A 1 J i0 + 3A 2 (J 50 + J51PKVV 
-3(J 31 + [AoJsi - AiJa + A 2 (3J 51 + 5J 52 )}U) A^ 

+3(B J 41 - BMq^vTuV - 3(B J 42 - B^) A^ g A ) + 6J 52 C ir^u x \(3A3) 

and the expressions for the Ai, Bi and d (known for a given equation of state) gives a 
set of field equations for the variables n, e, u", II, and q" which contains only three 
positive valued functions of n, e, namely Cn, C q , and C n . The C^'s are the collision integrals 
that depend on the microscopic interactions such as cross-sections. The primary transport 
equations depends on the C^'s and the local distribution function f(x,p). The relaxation 
times and length then follow from the product of the primary and the second order transport 
coefficients. The resulting relaxational transport equations [7] are presented in in detail for 
applications in [1]. 

In summary the 14-Field Theory of relativistic non-ideal fluids is concerned with the 
conservation of net charge and of the energy-momentum and the balance of fluxes. We 
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rewrite the conservation and balance equations by splitting them into spatial and temporal 
parts and by making right hand sides explicit. We have 



d^N" = 0, u v d\T u = 0, A»d x T uX = 0, (3.44) 

UvUxdyF"* = -4C U A 2 U, Afr^F"* = C q B iq ", d x F^ x = l -CX^ v) ■ (3.45) 
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IV. 2 nd ORDER ENTROPY 4-CURRENT IN KINETIC THEORY 

The 14-Field theory is restricted by the requirement of hyperbolicity and causality. These 
requirements can be deduced from the entropy principle. The kinetic expression for entropy, 
Eq. (2.27), can be written as 

S" = -J dwp^ip(f) , (4.1) 

where 

W) = [f(x, p) In (A^f(x, p)) - a^AoAix, p) In A(x, p)] 

= -a- 1 A lnA(x,p) + f(x,p)\n[A^ 1 f(x,p)/A] . (4.2) 

Expanding ip(f) around ip(f eq ) up to second order, i.e., 

w) = v>(/ eq ) + w eq x/ - n + \r(nu - rf + ... , (4.3) 



where 



= \n[A 1 f(x,p)/A(x,p)} , (4.4) 
r(f) = [A(x,p)A- 1 f(x,p)}- 1 , (4.5) 



are the first and second functional derivatives of tp(f) with respect to f(x,p), gives 

= -a-'AolnA^p) + [a(x) - (3^x)p^[f(x,p) - r(x,p)\ 

+i [r(x,p)A- 1 A^(x,p)}- 1 (f(x,p) - r(x,p)) 2 + ... . (4.6) 

Inserting the expression for ip(f), Eq. (4.6), into the expression for entropy flux, Eq. (4.1) 
yields 

s» = s^ + q ^- l -j dw p» nnif - r? , (4.7) 

where 

S£ = P (a, (3)(3» - aN^ + j3 x T% , (4.8) 
13 



is the equilibrium entropy and is obtained by inserting f eq (x,p), Eq (2.2), into Eq. (2.27). 

Using Eqs. (3.26)-(3.28) for e,e M , and in terms of the fluxes, one substitute f(x,p), 
Eq. (3.1), into Eq. (4.7) to get the non-equilibrium entropy 4-current, 

S» = + ^ - ±p U n (ftll 2 - [ft + w- 1 ] q X q x + (3 2 n»\ vX ) 

-(3 ( [a + w- 1 ] <f n - + w' 1 ] q x ^ x ) . (4.9) 

The coefficients a i: ft stand for [cf. [7]] 

-P41-P20 — -P31-P30 rzLirn 

qj _ ^31 1^52 — Ja\Ja2 ^ 

3 f 3 1 

^° = /9 t2 02 1 '^ 52 — p^[<^3i (1/31^30 — J41J20) + 1/41 (i74ii7io ~~ ^31^20)] f ? (4-12) 

A = ^2^ , (4-13) 

1 J 5 2 

2ft7 4 2 



& = i^lr • (4-14) 



These seccond order coefficients involves the J n k which can be obtained as differentiations 
of the X n k (or the equation of state p = p(a, ft) with respect to a and ft The transition to 
p = p(e,n) is effected by the relations 



dp 
oa 



(4.15) 



With the help of Gibbs' equation 

d(pft) = JV^da - T^dft , (4.16) 

the entropy production can be written immediately as in phenomenology. In kinetic theory 
we can derive it by substituting Eq. (3.1) into the expression for entropy, Eq. (4.7). Then 
we invoke the second law of thermodynamics, d^S^ > 0, and developing up to second order 
in the deviation function. This leads to an expression for Td^S^ which can be written 
as the sum of contribution Td^S^ linear in a deviation function <f>(x,p) and another one 
Td^S^ which is quadratic in <fi(x,p). Using Eqs. (3.26)-(3.28) for e, e M , and in terms of 
dissipative fluxes together with the help of the expression for the derivative of f eq (x,p) 

p^(lnA ^D = = -pPd, {^^} , (4.17) 
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we obtain the following contributions 

Td^s? = -n Vm ^ - ?" " D ^) + ^ v ^ UlJ) ' (4 ' 18) 

- Vm { [<*o + ^ - [a x + w" 1 ] ^tt""} . (4.19) 

The second law of thermodynamics requires the form of he collision term C[f] to be such 
that the entropy production is given by a positive-definite integral for any function f(x,p). 
That is, 

J dw p»d^(f) = -dpS" = j dw < . (4.20) 

Eq. (4.20) together with the first two moment equations, Eqs. (3.36) and (3.37), and 
Eq. (3.19) leads to 

< = - J dw C[f]y(x,p) = -e^Pi" (4.21) 
The entropy production can be written as 

/T 1 ^ = C~ X n 2 - \~ 1 q a q a + (2r ] )- 1 7r a ^ > , (4.22) 

with 

^vjiJMt, ^ KT = JL = ^, „ = JL = 10 ^ (4 . 23) 



where the Ca are the collision integrals which involves the cross-sections of various processes. 
For the calculation of the primary transport coefficients see, for example, [9]. For consis- 
tency the primary transport coefficients, the second order transport coefficients (hence the 
relaxation times and relaxation lengths) and the equation of state have to be determined 
from the same model or theory. 

From the kinetic theory approach the unknown phenomenological coefficients 
(, k, rj, a , cci, Po, Pi, and f3 2 can now be explicitly identified from the knowledge of the 
collision term and the equation of state, k stands for the thermal conductivity, and ( and 
r) stands for the bulk and shear viscous coefficients respectively. The transport coefficients 
involve complicated collision integrals. The relaxation coefficients f3 , Pi and p 2 makes the 
theory a causal one. The coefficients a and cci arise from coupling between viscous stress 
and heat flux. Knowledge of the second order coefficients allows one to write the primary 
transport coefficients in terms of the relaxation times. Such relaxation times depend on the 
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collision term in the Boltzmann transport equation, and their derivation is an extremely 
laborious task [9, 10]. For present purposes, it suffices to know that the kinetic theory yields 
the form of the second-order entropy 4-current and of the evolution equations of the fluxes 
and that it provides the explicit values of the relaxation coefficients in them. The relaxation 
times are related to the transport coefficients multiplied by the (3a and and the relaxation 
lengths are related to the transport coefficients multiplied by the «a, 

r n = CA), r q = kT(3^ r n = 2r^ 2 , (4.24) 
lu q = C a o, = kToo, l qn = kT«i, l nq = 2i]ai . (4.25) 

These are the relaxation times for the bulk pressure (rn), the heat flux (r q ), and the shear 
tensor (r n ) ; and the relaxation lengths for the coupling between heat flux and bulk pressure 
(Inq, lqn) an d between heat flux and shear tensor (l qn , l nq ). From these one can already study 
the general dependence of the ratios of the transport coefficients to the respective relaxation 
times or relaxation lengths, 

(4.26) 
(4.27) 

where \ qq stands for kT. The various time scales and length scales mention here are to 
be contrasted with appropriate dynamical time scales and the transverse and longitudinal 
length scales for a given nuclear collision. In high energy nuclear collisions the time scales 
over which the various constituents of the matter produced in these collisions equilibrate 
are important in understanding the nonequilibrium effects. Thus a knowledge of relevant 
nonequilibrium properties is essential for a complete description. In the collision dominated 
regime the mean collision times and the relaxation times associated with the changes in 
distribution functions provides the information about the trends towards global dynamics. 
These times are generally smaller than the times associated with the transport of momentum 
and energy. 

Once the entropy 4-vector is obtained from the kinetic theory one may then construct a 
generating function which may be a 4-vector [11] or a scalar [12]. Thus from the entropy 
4-vector one constructs/define a 4-vector $ = ^(e', e' x , e v \) 

W = S» + e'N^ - e' x T^ + e uX F^ x (4.28) 
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where e' = a + e and e' A = j3\ + e A . From the differential form of ^ 

= N n de ' _ T ^ de ' x + F^ A de^ A (4.29) 

and the entropy principle (the second law of thermodynamics) which is written as 

dpS" + e'd^N* - e' x d^ x + e uX (d^ x - P uX ) > (4.30) 

one then obtains the functions N^, ', and F^ A ^ as functions of e', e' A , e u \. That is, 

^ M = = FM" = - -g» x g a °, -e, u P^ > (4.31) 

de' de\ de uX 4 de a p 

Using the definition of Eq. (4.28), we can convert Eq. (4.29) into a differential form for 
the entropy 4-vector 

dS"" = -e'dJV" + e^dT^ - e^dF""* (4.32) 

In discussing the constraints imposed on the dynamical equations by the entropy prin- 
ciple, hyperbolicity and causality requirements we need to cast the system of dynamical 
equations in a more transparent form. We introduce Y% = {N^, Tjf, F^ uX) } representing the 
primary dynamical variables, P% = {0, 0, P v \} representing dissipation source tensor and 
X A = {V, e' A , e u \} representing the auxiliary dynamical variables. A = 1, • • • , 14. Then the 
dynamical equations Eqs. (3.36)-(3.38) with the restrictions Eqs. (4.31) can be written in 
the form 

d,Y% = P A , and n = ~^r (4-33) 
which can be combined to form a symmetric form 

-V** " 7^1^ Vm X b = P A (4.34) 



dX A dX B " D dX A 0X B 

For the system of equations Eqs. (4.34) to be hyperbolic (i.e., has well-posed initial value 
formulation) and causal (hyperbolic with no fluid signals propagating faster than light) we 
require that 

(4.35) 



dX A dX B 

be negative-definite for the physical states of the fluid. 
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The requirements of non-negative entropy production and of hyperbolicity imply the 
following restrictions on the coefficients Cn, C q , C n and f3 Q , (3±, (3 2 respectively. These read 



C n > 0, C q > 0, C n > 
A> > 0, A > 0, & > 



(4.36) 



(4.37) 




(4.38) 



The conditions Eq. (4.36) ensure the positivity of bulk viscosity, shear viscosity and heat 
conductivity. The conditions Eq. (4.37) ensure that the entropy density has its maximum in 
equilibrium and together with conditions Eq. (4.36) they also ensure the positive relaxation 
times for II, q^, and ir^ . The two inequalities in Eq. (4.38) are the stability conditions on 
compressibility and specific heat. The whole set of inequalities Eqs. (4.36)-(4.38) ensures 
finite speeds. Some of the consequences of the requirement of hyperbolicity impose the 
bounds on the values of II, g M , and terms of the independent variables of the equation 

of state e.g., p and e. The hyperbolicity requirement, that is the requirement that the field 
equations form a symmetric hyperbolic system ensures that the problem is well posed and 
that the characteristic speeds are real and finite. This requirement is equivalent to the 
requirement that the second differential 5 2 S tJ, (n, e, II, g M , ir^) of the entropy must be negative 
definite. 

V. THE EQUATION OF STATE PRESCRIPTION 

For proper description of the space-time evolution of the hot and dense nuclear matter, 
produced in high energy nuclear collisions, using fluid dynamics, one needs an equation of 
state to close the system of evolution equations. The pressure as a function of energy density 
and baryon density, p(e, n), is required for solving the fluid dynamical equations numerically. 
It is convenient to tabulate this function, since calculating the pressure in the hadron phase 
for given e and n requires a double root search to find T and fi. To facilitate fast and easy 
numerical computation one discretizes the e — n plane. Every other thermodynamic quantity 
can be calculated as a function of e and n. Of particular importance are temperature, 
baryochemical potential, entropy density, the relaxation and coupling coefficients. For the 
fluid dynamical calculations intermediate values of thermodynamic quantities are calculated 
by two-dimensional linear interpolation. 
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For this case study the nuclear matter is described by a a — 0-type model [13] for the 
hadronic matter phase and by the MIT bag model [17] for the quark-gluon plasma, with a 
first-order phase transition constructed via the Gibbs phase equilibrium conditions. This 
type of an equation of state was presented in [14, 15]. 

In the hadronic matter, the equation of state, i.e., the pressure p as a function of the 
independent thermodynamical variables temperature T and baryochemical potential p, for 
non-strange interacting nucleonic matter is defined by [16] 

Phad(T,fi) = p N (T,v;M*)+p N (T,-v;M*) + J2Pi( T ;mi) (5.1) 

i 

+ n V(n) - / V(n') dn' - p s S(p s ) + / S(p' s ) dp' s . 
Jo Jo 

Here, 

p N (T, v- M*) = J£(0*, z*) , p N (T, -v- M*) = z*) , (5.2) 

where 0* = (3v and z* = f3M*, is the pressure of an ideal gas of nucleons moving in the 
scalar potential S and the vector potential V. These potentials generate an effective nucleon 
mass 

M — ► M* = M — S(p s ) , (5.3) 

where M = 938 MeV is the nucleon mass in the vacuum, and also shift the one-particle 
energy levels 



E ± = Vp 2 + M 2 — > E* ± V(n) = \/p 2 + (M*) 2 ± V(n) . (5.4) 

The vector potential is conveniently absorbed in the effective chemical potential 

i/ = n-V(n), (5.5) 

giving rise to the interpretation of Eq. (5.2) as the pressure of an ideal gas of quasi-particles 
with mass M* and chemical potential v. Furthermore, 

Pi{T ]mi ) =F 21 (0, Zi ) , (5.6) 

where Zi = ftrrii, is the pressure of an ideal gas of mesons with degeneracy gi and mass mj. 
Only the pions (g n = 3, m n = 138 MeV) are considered since they are the lightest and thus 
most abundant mesons, n is the net baryon density, 



(rp \ dp 



= Jf o (0*,^)-Jf o (-0*,^), (5.7) 

T 
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and 



1 



e {E*-v)/T + I e (E*+v)/T + ]_ 



(5.8) 



is the scalar density of nucleons. 

Once the pressure is known, the entropy and energy density can be obtained from the 
thermodynamical relations, 



dp 

S= df 



e = Ts + pin-p . (5.9) 

The potentials V, S are specified as in Ref. [15] 

V{n) = C 2 v n-C 2 d n l '\ S(p s ) = C 2 p s , (5.10) 

where C\r = 238.08 GeV -2 , C| = 296.05 GeV" 2 , Cj = 0.183, are the parameters which 
which leads to the reasonable values for the effective mass and the compressibility of Mq = 
0.635 M, K = 300 MeV. T and \i emerge naturally from the double root search. Then 
using the results from Appendix A, i.e., Eqs. (All) and (A12) or Eqs. (A40) and (A41), one 
calculates the thermodynamic integrals X n k and J n \. for given e and n Then the aii(e, n) and 
fii(e,n) are determined. 

In the quark gluon plasma phase we employ the standard MIT bag model [17] equation 
of state for massless, non-interacting quarks and gluons, i.e., 

PQG p(T,ri=l±(0,(i>)-B , (5.11) 

where is given in Appendix C and B is the bag constant. Other thermodynamical 
quantities follow again from Eqs. (5.7) and (5.9). Note that p does not depend explicitly on 
n for this EoS, 

p= l -{Z±-AB) . (5.12) 

The value of the bag constant is taken to be B = (235 MeV) 4 which results in a phase tran- 
sition temperature of T c ~ 169 MeV at vanishing baryon density. Other thermodynamical 
quantities again follow from Eq. (5.9). The baryon chemical potential hb is related to fi q 
by fis = 3/i g and the net baryon charge density is related to the quark-gluon plasma 
baryon charge density by n B = rib/ 3. In the quark-gluon plasma the n = case is simple 
since we also have ix — 0, and one obtains from Eqs. (5.11) and (5.12) the simple formula for 
the temperature, T = [60e/(32 + 21 N^ir 2 } 1 / 4 . Nf is the number of quark flavours. In our 
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equation of state we consider Nf — 2. For finite //, T can be eliminated from the equation 

of e using the equation for n. This results in a sixth order equation in fj,, 

6N f -8 B 21A f - 58 o 9 (32 + 21A f )7r 2 9 

°=i2k^" — ^""^"-^Sir" • (5 ' 13) 

which has to be solved numerically. After that, T = [9n/Nffi— [i 2 /Qn 2 ] 1 / 2 (which follows from 
the equation for n). Once T and ji are known for given e and n, the thermodynamic integrals 
I n k and J7„fc are also calculated as functions of e and n using the results of Appendix C and 
thus the second-order coefficients aii(e,n) and f3i(e,n) are also calculated. 

In the mixed phase the quark gluon plasma phase equation of state Eq. (5.11) is matched 
to the hadronic phase equation of state Eq. (5.1) using the Gibbs' phase equilibrium condi- 
tions, 

Phad = PQGP, Thad = TqgP, ^had = f^QGP ■ (5-14) 

In the mixed phase, for given temperature T and chemical potential ji the values of energy 
and baryon density read 

e = ^Q £ q(T, A*) + (1 - \ Q )e H {T,n) , (5.15) 
n = \ Q n Q (T,n) + (1 - X Q )n H (T,fi) , (5.16) 

where \q is the fraction of volume the quark gluon plasma occupies in the mixed phase. 
Conversely, for given e, n these equations yield values for Xq, T, jj,. This is done numerically 
using the values of Eh,q{T, p), uh,q(T, /j,). Once T and \x are known the pressure follows 
from Gibb's phase equilibrium condition, Eq. (5.14). Similarly the other thermodynamic 
quantities can be calculated as function of e and n. Of particular importance are tempera- 
ture, baryo-chemical potential, entropy density and the rest of the thermodynamic integrals, 
X n k and J n ki from which the second-order coefficients aii(e,n) and f3i(e,n) are evaluated. 

In Figs. 1 and 2 we show the dependence of the relativistic thermodynamic integrals on 
this particular equation of state. The results are only for the hadronic phase without the 
phase transition into the quark gluon plasma. The results of the latter scenario will be 
explored in full detail somewhere else. 

VI. FREEZE-OUT AND PARTICLE SPECTRA 

At any space-time point of the many-particle (fluid dynamics) evolution a particle 
"freezes out" when its interactions with the rest of the system cease. This depends on 
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the mean free path of the particular particle. If the mean free path is small enough such 
that local thermodynamical equilibrium is established, fluid dynamics is valid. For larger 
mean free path non-equilibrium effects become increasingly important, and if the mean free 
path exceeds the system's dimension, the latter starts to decouple into free-streaming parti- 
cles. To describe the system's evolution in the latter two stages in principle requires kinetic 
theory. 

The region where the mean free path is of the order of the system size can be approximated 
by a hypersurface in space-time [18, 19]. When a fluid element crosses this hypersurface, 
particles contained in that element freeze out instantaneously. The freeze-out should be 
treated in a self-consistent way taking into account the non-equilibrium effects caused by 
particles leaving the system. 

We will consider isothermal freeze-out which corresponds to freeze-out at a constant fluid 
temperature Tf. One could also consider freeze-out at a constant center of mass time tf. In 
the later scenario the temperature along the hypersurface is no longer constant. 

In the calculations of particle spectra at freeze-out, the distribution of particles that have 
decoupled from the fluid is determined as follows. The Lorenz-invariant momentum-space 
distribution of particles crossing a hypersurface £ in Minkowski space is given by [19] 



where dS M is the normal vector on an infinitesimal element of the hypersurface E. d£ M is 
naturally chosen to point outwards with respect to the hotter interior of E since Eq. (6.1) 
is supposed to give the momentum distribution of particles decoupling from the fluid. The 
distribution function f(x,p) is given by Eq. (3.1). 

The calculations simplifies in the case of longitudinal boost invariance and transverse 
translation invariance [20]. In this case the on-shell phase space distribution function de- 
pends only on the reduced phase space variables (r, x, P±), where r 2 = t 2 — z 2 and \ = 1 l~Ui 
in terms of the (longitudinal) particle rapidity y = tanh" 1 \p z / 'p°) and the (longitudinal) fluid 
rapidity rj = tanh _1 (z/t). The collective flow velocity field, the momentum and the heat 
flux in this case are given by 




(6.1) 



w 



p 



x^/t = (cosh?7, 0j_, sinh rj) , 

(m_i_ cosh y, p± cos <f> p , p± sin P , m± sinh y) , 

gP = g(sinh?7, 0j_, coshr?) , 



(6.2) 
(6.3) 
(6.4) 
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where / M is a space-like 4- vector IH^ = — 1 and it is orthogonal to the 4- velocity l^u^ = 0, 4> p 
is the angle of particle's momentum around the z-axis and m± = \Jp\ + m 2 is the transverse 
mass which reduces to p± for the massless particles. The shear tensor is given by 



7T 



f 7i s sinh 2 r) 



\n s cosh 7] sinh 77 



2 








7i s cosh r) sinh 77^ 



2 u 

7T S cosh 2 ?7 J 



(6.5) 



where 7r s is the scalar shear pressure (to distinguish it from the constant n). In our particular 
case at freeze-out r = tj, E is conveniently parameterized by S M = (ry cosh 77, xj_, 7/ sinh?]) 
with dS' 1 = (r/ cosh 77, 0j_, —r/ sinh r])dr]dx_L. Therefore, d£ M p M = r/mj_ cosh xd7/dxj_ and 
j9^-u M = mj_ cosh x- 

The distribution function however now includes the corrections due to dissipative ef- 
fects. Now there is the equilibrium part and non-equilibrium part. With the expressions, 
Eqs. (3.26)-(3.28), for e(x), e u (x) and e^ v (x) in terms of the dissipative fluxes we can write 
the deviation function as 

<j>(x,p) = [Ao ~ A(p%) + [4(j9%) 2 - m 2 ] A 2 ] n+[/3 - Bi(p%)] (p%)+C p fl p u rr^+- . 

(6.6) 

In our present freeze-out prescription we have 



pfq^ = qm± s inh \ , 
V^f^nv = Tt s m\ sinh 2 x ~ \ n sP\ , 



(6.7) 
(6.8) 



and the deviation function becomes 



<j)(x,p) = [A - Ai(m± coshx) + [4(m_L cosh x) 2 - m 2 ]„4 2 ] II 
+ [Bo — Bi(m± coshx)] (771 j_ sinhx)? 



+C 



m\ sinh 2 \ ~ ^P± 



7T. S 



(6.9) 



In the freeze-out prescription of a system with longitudinal boost invariance and transverse 
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cylindrical symmetry [21] one frequently uses the integrals of the form 



p2n poo 

Xijmn = / d^sin>cos J V / dx sinh m x cosh™ x 
Jo Jo 



x c — r~A> > ( 6 - 10 ) 

g2x cosh x~ a ± cosip— q> q 1 x ' 



p2n /»oo 

Xifmn = dip sin* ip cos J ip dx sinh m \ cosh™ \ 
Jo Jo 



„zj_ cosh x~ai cos ip— <fr 

x 9 , (6.H) 

|gZ± cosh x~ a ± cos ip—<t> qJ 

where ip = — <p p is given in terms of the azimuthal angle of the surface element of the fluid <p> 
and the angle of the particle's momentum <p p around z, z± = r y±m±/T and a± = r y±v±p±/T. 
The powers in the trigonometric and hyperbolic functions comes from the product of the 
volume element parameterization in cylindrical coordinates with boost invariance and the 
distribution function f(x,p) = fip^u^, Il,p tl q ll ,p tl p v 7r ,J ' 1 ') as given by Eq. (3.1). Eq. (6.10) 
arises due to the first term of Eq. (3.1) and this gives the equilibrium (ideal fluid) contribution 
to the spectra while Eq. (6.11) arises from the second term of Eq. (3.1) and it gives non- 
equilibrium (non-ideal fluid) contribution to the spectra. Expanding the Fermi-Dirac and 
Bose-Einstein distribution functions in geometric series we can write Eqs. (6.10) and (6.11) 
as 

oo 

Xijmn 

fc=l 

/*2"7T pOO 

x d^sin^cos^e^ 08 ^ / d X sinh m x cosh™ X e~ kz± coshx , (6.12) 
Jo Jo 

oo 

xgL = EM' -1 *** 

k=i 

P'2tt poo 

x dV>sin l >cos J >e fea±cos M dxsinh m xcosh"xe- fczxCOshx , (6.13) 
Jo Jo 



where the upper sign is for fermions and the bottom one is for bosons. Using the integral 
representation of the modified Bessel functions of the first kind I n (x) and of the second kind 
K n {x) we can write 



X 



(i) 

ijmn 

k=l 

oo 

(2) 



J2 (T)"' 1 e k %(ka ± )k mn (kz ± ) , (6.14) 



xSmn = ^2 (t)^ -1 ke k(f> iij(ka±)k mn (kz±) , (6.15) 



k=i 
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where 



[(1/2H / \(l/2)d]\ 

% = J2 y (2b + 2l-l)\\y- b - l 2nl b+l+h (y) , (6.16) 
[(1/2)4 / r fl/2W1 \ 

L= ( / 1 (2fc + 2Z-l)!!y- 6 -%+^(y) , (6-17) 

and [p] stands for the largest integer not exceeding p. In iij, (d — j, b — i/2, y — ka±) while 
in k mn , (d — n, b — to/2, y = kz±). Furthermore, h is zero or one depending on whether 
d is even or odd respectively. Note also that for odd powers in sin^ the vanishes and 
similarly for odd powers of sinh \ the k mn vanishes. 

In the case of pure longitudinal expansion i — j — and a± = 0. Then iij(ka±) = 
ioo (0) = 2n and Xijmn reduces to 

oo 

X&L = 2vr (T)^ 1 e^AvnnC^x) , (6.18) 
k=l 

oo 

X$L = ^J2 to*"* ke k *k mn (kz±) • (6.19) 
fe=i 

The inclusive particle distribution and transverse energy at freeze-out are given by 
dN(r f ) 



dyd 2 p ± 
dEi 



J Tfdi]d 2 x±m± cosh x f(x,p) , (6.20) 

/ d 2 p±m ± dN }^ = ! r / d? 7 d 2 x ± I d 2 P± micoshx/(x,p) . (6.21) 
dy J dyd p± J J 

Using Eq. (3.1), for the distribution function f(x,p), the particle spectra can be written as 

^5f^ = [ T f dr ) d 2 x ± m ± cosh X {n(x,p) + r(x,p)A e *(x,p)<j>(x,p)} , (6.22) 
dyd p ± J 

where 

™ = - ^* p{m± co S h x /r } -i - < 6 - 23 > 

The first term in the curly brackets of Eq. (6.22) is the equilibrium contribution to the 
spectra while the second term is the non-equilibrium contribution to the spectra. Using the 
expression for <f>(x,p) and the integral representation Xijmn the particle distribution becomes 



dN(r f ) 1 



, - A -R 2 ± T f x \ xSi + (A - toM 2 )Xo 2 ooi - Aim±x<wo2 + ^mlxoow 
dyd p ± 2 



2 (2) 1 2 (2) 



n 

tt s !> (6.24) 
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where the first term which involves the xfjinn ls the equilibrium spectrum and the rest of 



the terms that involves Xi/mn represent the non-equilibrium contribution to the spectrum, 
i.e, they represent the corrections to the particle spectra due to dissipation. For dissipative 
corrections of particle spectra due to shear viscosity only we find, using Eqs. (6.10) and 
(6.11), 



dN(r f ) 



AoirR ± Tfm,j_- 



7T= 



a 



K 3 (z ± ) 



1 



(6.25) 



dyd 2 P± ^' J -^Ae + p 

In the case of a pure (l+l)-dimensional expansion in planar geometry, the transverse 

dimension enters only as a (constant) transverse area factor A. The parametric integration 

over the hypersurface £ in Eq. (6.1) is performed by distinguishing space-like parts S s (with 

z as integration variable) and time-like parts E t (with t as integration variable). This yields 

the total momentum-space distribution 

dN 



dyp±dp A 



= 2nAm^ 



dz [coshy — smhy(dt/dz)z\f(x,p) 



+ 



/ dt [cosh y(dz/dt)z — sinhy]/(x, p) 



(6.26) 



where (dt/dz)^ is the (local) slope of the space-like hypersurface element. In the expression 
for the distribution function Eq. (3.1) the equilibrium distribution f eq is given by Eq. (6.23) 
while <p(x,p) is given by Eq. (6.9). Note that in pure (l+l)-dimensional expansion we have 
only one independent component of heat flux which we take to be q z = 7Q 2 = 75 and one 



independent component of shear stress tensor which we take to be n zz = ■y 2 r zz 



7 2 7T S . 



We simply write q for the local rest frame independent component of heat flux Q z and n s 
for the local rest frame independent component of shear stress tensor r zz . From [1] we 
know that the other nonvanishing component of heat flux can be written as q° = q z v z and 
the other nonvanishing components of the shear stress tensor are ir 00 



Tl zz vl 



IT 



Oz 



7T 



zO 



7r s 7 2 f 2 and n xx = tt 



yy 



-7T s /2. 



We now apply Eq. (6.26) to the isothermal freeze-out scenario. Along the isotherm, 
T = Tf = constant and only i] depends on position or time, respectively. The rapidity 
distribution is obtained by integrating Eq. (6.26) over transverse momentum. For massless 
particles, 

= 2nA 

; , t coshy(dz/d^-sinhy^ ^ 



dy 



^ coshy — sinhy(dt/dz)-£, 
cosh 3 x(z) 



cosh x(t) 
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where 



T = — [Jio + [J 10 Ao - J 20 A 1 + 4J- 3 o]n + [B J 20 - B^qtemhxiz) 

47T 



(6.28) 



2 cosh 2 x(z) 

and the T nk and the J nk are given in Appendix C for single particle species. Analogously, 
the transverse momentum distribution is obtained integrating over y, 

dN 



where 



and 



P±dp A 



2tt p ± 



/oo 
dy [cosh y — sinh y(dt/dz)s] J 7 
-oo 

oo 



+ 



dt dy [cosh y(dz/dt)-z — sinhyjjF 



T = A 



exp{(p ± /T f ) coshx{z)} - a 



exp{(p ± /T f )cosh X (z)} 
°[eM(P±/Tf)cosh x (z)}-a]^ P) 



(6.29) 



(6.30) 



0(x, p) = [A - Aip ± cosh x(z) + 4p ± cosh 2 x(z)A 2 ] II 

+ [B — Bip± cosh x(z)} p± sinh X (z) q 

1 



+C 



p\ sinh 2 x(z) - -p\ 



(6.31) 



In the isochronous freeze-out scenario the hypersurface X has no time-like part and the 
second term in Eq. (6.26) vanishes. Moreover, also the second term in the numerator of the 
remaining term vanishes due to dt/dz = for this particular hypersurface. The temperature 
along the hypersurface, however, is no longer constant. Thus the rapidity distribution 
becomes 

dN n 2 , f L+tf , coshy r 

— = 2n 2 A / dz y T , 6.32 

dy J- L -t f cosh 3 x(z) 

with T{z) given by Eq. (6.28) with Tf now replaced by T(z). The transverse momentum 

distribution is 

dN rL+tf POO 

— = 2nA p± dz dy coshy T . (6.33) 

P±dp± J-L-t f J-oo 

where T is given by Eq. (6.30) with the temperature Tf now replaced by T(z). The results 
of this section generalizes that of Ref. [22] by including non-equilibrium (dissipative) effects 
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in the freeze-out prescription. The most recent applicatio of freeze-out prescription in the 
presence of shear viscosity corrections has proven to be excellent in describing the spectra 
of particles produced at RHIC [23]. 

VII. SUMMARY AND CONCLUSIONS 

We have presented a microscopic interpretation of dissipative fluid dynamics by means 
of the 14-moment method (or 14-Field Theory: an effective kinetic theory for relativistic 
non-equilibrium fluid dynamics). From the relativistic kinetic theory using Grad's 14 mo- 
ment method one obtains relativistic causal fluid dynamics which is a field theory of the 14 
fields of net charge density-particle flux and stress-energy-momentum. The field equations 
are based on the conservation laws of net conserved charge, energy-momentum and on a 
balance of fluxes. The equations satisfies the principle of relativity, entropy principle and 
the requirement of hyperbolicity (for causality). The resulting field equations contain only 
bulk viscosity, shear viscosity and heat conductivity as unknown functions. All other coef- 
ficients may be calculated from the equilibrium equations of state. The interface between 
the microscopic and microscopic in the causal dissipative fluid dynamics is effected via the 
standard transport coefficients. The comparison between the microscopic and the macro- 
scopic descriptions provides more information on the transport properties that govern the 
relaxation of various dissipative processes. 

From the kinetic theory approach the unknown phenomenological coefficients 
C, k, T), a , a±, j3 , Pi, and f3 2 can now be explicitly identified from the knowledge of the 
collision term and the equation of state. We conclude that the 14-Field Theory of viscous, 
heat-conducting fluids is quite explicit - provided we are given the thermal equation of state 
p = p(a, (5) -except for he coefficients Ca- These coefficients must be measured or determined 
from the underlying theory/model which gave the equation of state. 

We have also presented the relativistic thermodynamic integrals in a more general for- 
malism and also for different special cases which are relevant for applications to nuclear 
collisions applications. These integrals are needed for calculating the second-order relax- 
ation coefficients. The results of the later will be presented in detail somewhere. 
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APPENDIX A: RELATIVISTIC THERMODYNAMIC INTEGRALS 



In relativistic kinetic theory one frequently encounters the moments of the distribution 
function describing local and global equilibrium. The nth moment is defined by 

( x ) = J dw p^...p^f(x,p) , (Al) 

and the corresponding change of the nth moment due to variation of f(x,p) is 

Sl^-^ix) = J dw p" 1 ...p fXn SyA(x,p)f(x,p) (A2) 

where y = ln(f(x,p)/A(x,p)). The variations of moments involves the auxiliary moments 

JK-> x »{x) = J dw p> Xl ...p t * n A(x,p)f(x,p) . (A3) 

These moments can be expanded in terms of symmetrized tensor product of and the 
metric tensor g^ u . We define for 2k < n, the set of rank n tensors 

/\(?k u n-2k) = 2 \k\(n — 2k)\ a^ 1 ^ 2 ^M2fc-i^2fe u V2k+i u ^n (A4) 

permutations 

where A^ v is the projector g^ v — u^u u . The sum, which runs over all distinct permutations 
of the tensor indices, has been divided by the total number of those permutations. The 
number k can take all integer values between zero and the largest integer not exceeding 
n/2. The latter value will be denoted as [n/2]. These tensors possesses the orthogonality 
property 

a (2k n 2k) a (2k + l)\(n - 2k)\ . ... 

A (2k u n 2k) = V J_\ J_S ql . (A5) 

71/. 

which arises because of all different terms, differing only by the trivial permutations like 
-> A"", A^ a A u/3 -> A u/3 A^ a , u»u v -> u u u^, only 2 k k\(2k)\ terms survived. If these 
trivial repetitions are lumped together the expansion contains 

a, lk = I " ] (2k - 1)!! (A6) 

terms. The double factorial notation stands for (2k — 1)!! = 1 - 3-5 • • • (2k — 1) for k — 1, 2, • • • 
and the recursion formula (2k — 1)!! = (2k + l)!!/(2/c + 1) extends the definition to negative 
integer arguments. The moment integrals are then expanded in the form 

[n/2] [n/2] 

j" 1 -"- = Yl a ^ z ^ A(2k un ' 2k) > = E A(2fe M "~ 2fe) • ( A7 ) 

k=0 k=0 
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The scalar coefficients T n k and J nk , which depend on the parameters a and (5 are found by 
contracting both sides of Eqs. (Al) and (A3) with a tensor of the form (A4) and using the 
orthogonality property (A5). The results of such contractions are 

A f°° 

Z nk (a,(3) = (2fc+ ° 1)! , J dw[p a u a - P y a ] k (p^r- 2k 

X e /3p"u„-a _ ~ ' ( A8 ) 

A f°° 

Jnk(a,P) = ( 2 fc + °i)!! J dw iP au « ~ P a P*\ k (P^v) 



x ? . (A9) 



fi-p v u v ~ot 

\£w" u »- a - a\ 

The above integrals are invariant scalars, therefore they can be evaluated in any frame. 
We will evaluate them in the local rest frame. In this frame u 11 = (1,0,0,0), so that 



p^Pn = p 2 = m 2 and p^u^ = p° = a/p 2 + m 2 . We introduce spherical polar coordinates 
d 3 p = p 2 dpdQ. Integrating over dfi and introducing new variables 

* = f, z = j, (MO) 

the above thermodynamic integrals, Eqs. (A8) and (A9) may be written as 

^A 2 r 
(2k + l)\\ J 

X e V*T^-* ' (AU) 



DC 



X nk {^z) = J nA " u T n+2 I dxx 2 ( fc+1 V + * 2 ) (n_2fc_1)/2 



n-2fc-l)/2 

i ^ y 

Vx 2 +2 2 — (f) 



X 2 ( A12 ) 

1 n-2k 

~P n—l,k—l ^ ^ -Ln—ljk ■ 

The second line of Eq. (A12) is obtained by partial integration. 

Derivatives of T n k((p,z) and J^nk{4>, z ), Eqs. (All) and (A12) can be expressed in terms 
of the Jnk{4>,z): 

dXnfc = Jnkd(p - m^^n+i^dz 

= Jnfcda - J n+ i,kd(3 , (A13) 
^d,7 nfc = m [J n -i,k-i + {n- 2k)J n - lik ]d<i> - [J n ,k-i + (n + 1 - 2k)J njk ]dz , 
P dJnk = [Jn-i,k-i + (n- 2k)J n ^ k ]da - \J n ^ + (n + 1 - 2A;) t 7„ 1 k]d/3 . (A14) 
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For the equation of state p = p(e, n) we need the derivative of p, e and n with respect to (0, 
z) or with respect to (a, (3). From Eq. (A13) we have 



dn = dJio = Jioda - J 20 dp, de = dJ 20 = J2oda - J 30 d(3 , 
dp = dJ 2 i = J 2 i da - J 3 id(3, dJ 3 i = Juda - Jndf3 . 



(A15) 
(A16) 



From Eqs. (A15) and (A16) the differentials of a(e ) n), f3(s, n), and p(e ) n) are written as 
da = {J 30 dn - J 20 de} , d/3 = — {J 20 dn - J w de} , (A17) 



dp 



1 



D 



20 



{(J21J30 - J31 J2o)dn - (J21J20 ~ JsiJio)de} . 



(A18) 



The fundamental thermodynamic relation together with the first law of thermodynamics 
can be expressed as 



nTd(s/n) = de- hdn = {J 20 - hj 10 )da - ( J m - hj 20 )d[3 . 



(A19) 



From Eqs. (A15), (A16) and (A19), the differentials of a(n, s/n) and (3(n, s/n) can be written 

as 



da = 



£20 
1 



{(J30 - hJ 2Q )dn - J 20 nTd(s/n)} , 



d(3 = — {(J 20 - hj 10 )dn - J 10 nTd(s/n)} , 
U20 



(A20) 
(A21) 



and the differentials of a (e, s/n) and (3(e, s/n) as 

1 



da = 
d/3 = 



hD 20 
1 



hD 



{(Jm-hJ 20 )de-J 30 nTd(s/n)} , 
{(J20 - hj 10 )de - J 20 nTd(s/n)} . 



(A22) 
(A23) 
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Thus from Eqs. (A16), (A20) and (A21) we obtain the expression used in bulk pressure 
relaxation coefficients, namely 



d In X 
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(9 In J10 



din J 4 



12 



s/n C> mJ 21 



J2lJ_ 

J42 D. 



20 



(A24) 



From the definitions of specific heats per net conserved charge 



C p = -{3 



d(s/n) 

df3 



,C V = -P 



d(s/n) 

df3 



, r — Cp/Cy , 



(A25) 
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one finds, with the help of Eqs. (A15), (A16) and (A19), 



c p - c v = p 2 (hj 10 - j 20 ) 2 /M 



10 J 



r- 1 



(hj w - JwY/D, 
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From the definitions of the speed of sound and the compressibilities 



2 dp 

c °=dS 



s/n 



dp 
Ik 



n dp 



e + p dn 



a, = 



1 dn 
n dp 



Kt = 



s/n 



1 dn 
n dp 



we find, with the help of (A15), (A16) and (A19), 



a. 



10 



10 



(A26) 
(A27) 

(A28) 
(A29) 



s (3{e + p)J 10 ' s n 2 T ' w 
Thus for a given equation of state, i.e. pressure as a function of two independent state 
variables, the other thermodynamic quantities can be obtained as partial derivatives of the 
pressure with respect to a and (3 or to <fi and z. This means also that the second order 
coefficients are also determined from the equation of state. 

From the net charge and energy conservation equations (in the particle or Eckart frame) 



n + n6 = , 

£ + (e + p + n)# + (v M <f + vr^ v^O = ° , 



(A30) 
(A31) 



with the standard notations A = u^d n A and 6 = <9 M w M , we can solve for a and (3 by 
simultaneously solving 



(A32) 
(A33) 

(A34) 

D - ' D ■ (A35) 

The l nk and J nk integrals, Eqs. (All) and (A12), can be written in terms of the more 
familiar functions ]C n (<f>,z), C n +i(<f>,z) defined, for n > 0, by (cf. [7]) 



J w a-J 2 o$ + IioO = , 

Jn<* - JaoP + (^20 + X 21 + 11)0 + (V M ^ + ^ v Vu u») = , 

where we have used Eq. (A15). We then find 

= J20X10 - Jio(X 2 q + X21 + n) g | Jiois/^ + tt^ y, ty) 



20 



20 



a = 



30X 10 ~ J 20 (X20 + X21 + II) e | JgojV^ + 7T^ <) 



K. n ((f),z) 
£ n+1 ((f),z) 



1 



(2ra- l)!!* n 
1 1 



2\-l/2 



dxx 2n (x 2 + z 2 ) 

e ^ x 2 +z 2_ ( f, _ a 

dx x 2n 



(2n-l)\\z n + 1 J Q e Vx^-<p_ a 
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(A36) 
(A37) 



Partial integration of these functions leads to the following relations 

^/C n = C n , -^C n+1 = -z n ^- {z- n lC n } = /C n _! + (2n/z)lC n . (A38) 

In the special case of a Boltzmann gas (a = 0), these functions becomes the Bessel functions 
of the second kind up to a factor of the chemical phase, 

£ n (<f>, z) = /C„(0, z) = K n {z) , (A39) 

where K n (z) are the modified Bessel functions of the second kind. 

By making use of the binomial expansion the integrals I n k and J7„fc, Eqs. (All) and (A12), 
can be expressed in terms of the K, n and C n+ i, Eqs. (A36) and (A37). For n = 0, 1, and 
k < this yields 

-n—k ( 

n V-^ ci.i.i', K-k+i+i for even n , 

X nk = 4rcA m n+2 ^ a nkl z^ k+l+ ^ I k+l+1 (A40) 

i=o ( Ck+i+2 for odd n . 

and differentiation of these results with respect to 0, (cf. Eqs. (A13) and (A38)), then gives 

, o 2 -*-~y n , ; , i\ I £>k+i+i for even n 

J nt = 4^ om »« J2 a na (A41) 

/=o ( >Ck+i + ^— z — -/Cfc+z+i tor odd n . 

The numerical coefficients are 

(2k + 2l + l)\\ (\n-k\ 

°»« = (M + ( 2 , J ' (A42) 

Note that in order to calculate the thermodynamic properties of ultra-relativistic (massless) 
particles using Eqs. (A40) and (A41) one has to take the limit m — > (hence z — > 0) of 
Eqs. (A36) and (A37). This is however not necessary if one uses the integral representation 
Eqs. (All) and (A12) since the case for m = or z = is included. The integrals T nk and 
J n k, Eqs. (A40) and (A41), satisfy the recurrence relations 

In +2 ,fc = rn 2 I nfe + (2k + 3)I„ +2ife+1 (2k < n) , (A43) 

where I; m stands for T[ m or Ji m . This follows from the contractions of Eq. (A7). 
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APPENDIX B: APPROXIMATION TO THERMODYNAMIC INTEGRALS 



1. Nondegenerate gas 

If cj) < z and z ^ 0, the integrals l nk and J nk , Eqs. (All) and (A12), may be evaluated 
[24] by making the substitution x = zsinhx expressing the functions (e zcoshx ± 1) _1 as a 
geometric series in e zcoshx and integrating term by term. The integrals can be written in 
terms of the Bessel functions with the help of 



1 oo 

= B^- 1 ^^ , (Bl) 

k=l 

oo 

= jy^ke-^ , (B2) 

where the upper sign is for the fermions and the bottom one is for bosons. The integral 
representation of the Bessel functions of second kind can be written as 

,oo KV2)«fl / \(l/2)d] \ 

/ dxsinh 2b X cosh d X e^ coshx = E (26 + 2r - l)\\y~ b ~ r K b+r+h (y) , (B3) 

^° r=0 y r J 

where h is zero or one depending on whether d is even or odd respectively. Then the integrals 
l n k and Jnk-, Eqs. (All) and (A12), can be written as 

. a oo [(l/2)d] 

= M ^ W l7 Q r 9 7 P E(Tl) Z ~ 1 [^ + (-l) n e- i 1 E 6 ^ y-^Hr+fcCy) (B4) 

^ Z=l r=0 

4 j oo [(l/2)d] 

where we have abbreviated 



Jnkr 




(26 + 2r-l)H, b = k + l, d = n - 2k , y = Iz . (B6) 



r 



The superscript notation ± on the X n k and J n ^ indicates that the integrals are evaluated 
taking into account pair production. The factor 7 p , which is 1 or 1/2 depending on whether 
the particle is charged or neutral respectively, takes into account the contribution of neutral 
particles in calculating thermodynamic properties of a gas. 7q is the quantum number for 
the conserved charge and i q is 1 or depending on whether n is even or odd respectively 
while j q is or 1 depending on whether n is even or odd respectively. 

34 



2. Extremely degenerate Fermi gas 



This case is relevant to the study of the low temperature region in high energy nuclear 
collisions and to the study of thermodynamic properties of neutron and quark stars. If 
(j)^$> z, the integrals K, n and £ n+ i, Eqs. (A36) and (A37), may be expressed in terms of the 
Chandrasekhar-Sommerfeld asymptotic expansion. 

M<M ~ 7 ^77- r X 2n (x 2 -Z 2 )- 1 / 2 dx + 2-- l -—-(<f- Z ^-^ ( j ) 

nVr ' ; (2n-l)\\z n J K ' 12(2n-3)!!z" vr ; r 

+2—- — rr7~~( ( t )2 - ^) (2n ~ 5)/2 0{3 + (2n - 5)(0 2 - z 2 )"V 2 } 

720 (2n - 5)!! z nKr J ri v JKr J r/ 

+2 31?r6 - — (d) 2 - z 2 ) (2n ~ 7)/2 6 

+ 30240 (2n- 7)!!^ 10 ^ j 

x {15 + 10(2n - 7)(0 2 - z 2 )- 1 ^ 2 + (2n - 7)(2n - 9)(0 2 - z 2 )' 2 ^} (B7) 

£ »«(^)~(^p^ 2 - 2 ) (2 * +,)/2 

+2 nr> 1 T)!!^^ 2 - z2 ) <2 "" 1,/2 ( 1 + < 2 " - Di* 2 - ^rV] 

+ 72O(2n-3)!!z"+ li0 j 

x {3 + 6(2n - 3)(0 2 - z 2 )' 1 ^ 2 + (2n - 3)(2n - 5)(<f - ^T 2 4 } 

+2 31?r6 - Lf^ _ 2x(2n-5)/2 

30240 (2n-5)!!^+ llV ; 
x {15 + 45(2n - 5)(0 2 - z 2 )- 1 ^ 2 + 15(2n - 5)(2n - 7)(0 2 - z 2 )" 2 4 

+ (2n - 5)(2n - 7)(2n - 9)(0 2 - ^ 2 )~ 3 6 } , (B8) 

with xq = (p. Changing the variable x = z sinh x one can evaluate the integral 

i rxo rxo 
R n = - x 2n (x 2 - z 2 )' 1 ' 2 ^ = z 2n / sinh 2 " X d X , (B9) 
z Jo Jo 
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with the help of Eq. (1.412.2) of Ref. [25] and the properties of the hyperbolic functions. 
Here we list the first few expressions for the K. n 

- 2 7vr 4 0(20 2 + 3z 2 ) 



MM-cosh-W*) T (02 _ z2)3/2 m W _ , 2)7/2 
3l7r 6 0(80 4 + 4O0V + 15^ 4 ) 

(BIO) 



1008 (0 2 - z 2 ) 11 / 2 

7T 2 77T 4 



*MM - -r 2cosh ' 1 ^) + > 2 -^ 1/2 + y (02 -%V 2 + 12O(0 2 -, 2 ) 5 / 2 
3l7r 6 0^ 2 (40 4 + 3z 2 ) 
+ 1008 (0 2 -^ 2 ) 9 / 2 ' 1 J 

^ 2 /C 2 (0, *) ~ ±* 4 cosh" 1 ^/*) + ^0(0 2 - z 2 ) V 2 (20 2 - 5^ 2 ) + ^0(0 2 - z 2 ) 1 / 2 

77T 4 0(20 2 -3^ 2 ) 317T 6 0z 4 

360 (0 2 - 2 2 ) 3 / 2 1008 (0 2 - z 2 y/ 2 ' 1 j 

^ 3 /C 3 (0, z) ~ -1/ cosh" 1 ^/*) + ^0(0 2 - z 2 )V 2 (80 4 - 260V + 3z 4 ) + ^0(0 2 - z 2 f' 2 
7tt 4 0(40 2 - 3z 2 ) 31?r 6 0(80 4 - 200V + 15z 4 ) 

+ 360 (02-^)1/2 +1512Q (02 _ ^2)5/2 ' ( ) 

and for the £ n+ i 

,£ l( 0,,)~(0 2 -, 2 )V 2 «* ^ + 



6 (0 2 - 2 2 ) 3 / 2 120 (0 2 - z 2 ) 7 / 2 
31tt 6 ,2 2 (80 4 + 120V + ,2 4 ) 



(B14) 



336 (02 -^2)11/2 

2/. , , x 1/ ,2 2^/2 ^ (20 2 -^ 2 ) 77T 4 2 4 

, 2 £ 2 (0,,) ~ g^-^ + y^^ + ^T^^ 
31 ^ 4 (60 2 + , 2 ) 

+ 1008 (0 2 -z 2 ) 9 / 2 ' 1 J 



317T 6 



6 



(B16) 



1008 (0 2 - z 2 ) 7 /2 ' 

that are needed for the calculation of the X nk and J nk used in the text. 

For a massless Fermi gas the same calculation can be performed for the X nk and J nk 
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directly using Eqs. (All) and (A12)) to give 



-nk 



AnA 
{2k + 1)!! 

+2(n + l) 



rpn+2 



Ln+2 



7T 



n + 2 

2 77r 4 
+ 



n ~ 2 317T 6 n " 4 

+ 777^777^ 777 + 



J, 



nk 



12 n\ 720 (n - 2) ! 30240 (n - 4) ! 

171+1 



(2fc+ 1)!! 



T n+ > + 1) 



n + 1 



-2 ^n-l ?7r 4 3^6 0n-5 

_l_2r7 1 I 1 I 1 I l 

" ' 1 12(n-l)! 720 (n - 3)! 30240 (n - 5)! 



(B17) 



(B18) 



3. Transition region 



The representation of the relativistic thermodynamic integrals as series of modified Bessel 
functions is useful for computation in the non-degenerate case and the representation in 
the degenerate region gets better for much greater than z. A method suitable for the 
transition region <fi ~ z is needed for computing the values of the integrals in this region and 
for numerical checking of known results. 

The method adopted for the transition region is based on making the substitution x = 
t(2z + t 2 ) l/2 in the integrals. The integrals Eqs. (All) and (A12) becomes 



nk 



(2£; + l)H 



rpn+2 2 



Jo 



f2(k+l)( z + t 2)n-2k( 2z + t 2)k+l/2 dt 



ot 2 +z-q 



J, 



nk 



4ttA) +2 f°° t 2(fc+1) (z + t 2 ) n ~ 2k {2z + t 2 ) k+l / 2 dt 



(2k + l)\r ~J [ e * 2 +^-0_ a ]2 
while the function JC n and C n+ i, Eqs. (A36) and (A37) are written as 



(2n 



£n+l — 



(2n 



i i r 
i i r 



t 2n (2z + t 2 )"- 1 / 2 cit 

e* 2+z -^ - a 
t 2n (z + t 2 )(2z + t 2 ) n - 1 / 2 dt 

pt 2 +z~S _ n 



(B19) 
(B20) 

(B21) 
(B22) 



APPENDIX C: ULTRARELATIVISTIC THERMODYNAMIC INTEGRALS 



1. Single particle densities for massless particles 



We want to study the thermodynamic properties of an ultrarelativistic gas. Ultrarela- 
tivistic particles are characterized by vanishing rest mass. We start by looking at the single 
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particle densities. For ultra- relativistic or massless particles, z = m/T = 0, the moment 
integrals, Eqs. (All) and (A12) take the simple forms 

4vrA _ h 



1, 



nk 



J, 



nk 



{2k + 

4ttA 
{2k + 1)!! 



^0 rpn+2 f 

1)!! Jo 



OO 1 

71+1 _ 



dx x r 



T 



n+2 



dxx 



71+1 



D x-4> 
[e x -<t> - a] 



(CI) 
(C2) 



For the fermions, when the equilibrium quantities are characterized by ~ 0, we can expand 
the the integrals in powers of 0. We have 
47rA n 



7- 

-L-nk 



{2k + 

x 



^0 rpn+2 f 

1)!! Jo 



1 e 2 
+ 



dxx 



n+l 



+ ^0 2 (n + 1)- 



+ 



(C3) 



Jr 



AttA 



nk 



T-rn+2 



{2k + 1)! 
1 

x 



dx x n 



^ J.2 X 

-<p n- 



n—l „x- 



_e x ~ ( f+l ' T [ e x-<t> + i] 2 2 x n [ e z-4> + i] 2 
With the help of (cf. [26]) 



+ 



Jo 



/oo i 
dxx^ 1 -^ = (l-2-(" +1 ))r(n + 2)C(n + 2), 



dxx 



n+l 



= (l-2"")r(n + 2)C(n + l) , 



(C4) 

(C5) 
(C6) 



(e* + l) 2 

where T(x) is the Gamma function and ({x) is the Riemann's zeta function, (Note that the 
C function with even argument is analytically known), we get, for a gas of quarks, 
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with qq the quark degeneracy. For the antiquarks we replace by —0. For gluons, = 0, 
and with the help of (cf. [26]) 
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- = T(n + 2)C(n + 2) , 
T(n + 2)C(n+l) , 
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we have 
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with g G the gluon degeneracy. 

As an example, the number density, energy density and pressure of quarks are given by 
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and for the antiquarks one replaces by — 0. Similarly, the number density, energy density 
and pressure of gluons are calculated using the boson integrals to give 
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ir' ~ 6U ~- 3 y 3 

We can also calculate the relaxation and coupling coefficients for the quarks and gluons. 
For = one can immediately check that for massless particles the quantity Q given by 
Eq. (3.33) that comes in the relaxation and coupling coefficients for bulk pressure vanishes 
identically. Thus those coefficients then diverges. That is, Q ==^ , «o =^ 00 , (3q 00. 
This have the consequence that the bulk viscosity itself has to go to zero. For the shear 
tensor relaxation coefficients, we have for the quarks and gluons respectively 
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4C(5)C(5)P 

For a classical Boltzmann gas in the ultra-relativistic limit (z<l) one gets the following 
expression for the shear relaxation coefficient 



(C19) 



2. A system of massless particles and antiparticles 



We now consider a gas of particles and antiparticles. For ultrarelativistic Bose gas the 
chemical potential is zero and the contribution of bosons to the thermodynamic properties of 
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the gas is given by Eqs. (Cll) (with n even, since for odd values of n there are cancellations, 
e.g., net baryon number vanishes) and (C12) (with n odd, since for even values of n there are 
cancellations). The same argument goes for the ultrarelativistic Fermi gas with vanishing 
chemical potential. But here we are considering a gas at finite chemical potential. For 
ultrarelativistic or massless particles, z = m/T = 0, the moment integrals, Eqs. (All) and 
(A12) takes the simple forms 
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These integrals can be evaluated by analytical means. We substitute y~ = x — in the 
first term and y + = x + in the second term and then write the integrals so that we can 
integrate from to oo. Then we have 
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For Jnfc we first perform partial integration. The first two integrals can be directly combined 
and the last two after the substitution y~ = —y + . We then have 
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In the last two integrals we noted that (e x + 1) 1 + (e x + 1) 1 = 1 and we have substituted 
z = x + (p. The two power functions in the numerator are then expanded according to the 
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binomial expansion. All terms of even j drops out: 
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The above integrals can be brought to their final analytic form by making use of 
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Thus we finally have 
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For the total particle number densities or single particle species densities separately one can 
no longer use the above formulation. One then uses the results from the previous section. 

To this end we give here the results of the ultrarelativistic thermodynamic integrals 
and J^ k for QGP. We present only those that are used in the text. We treat quarks and 
gluons as ultrarelativistic Fermi or Bose gases, respectively. Then from Eqs. (C29) and 
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(C30) we have , for a gas of quarks and gluons, 
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where go = N S (N^ — 1) and gq = N s N c Nf, (with N s the number of spin projections, N c 
the number of color charges and Nf the number of quark flavours) are the quark and gluon 
degeneracies respectively. Note that the net baryon charge is n = Xf Q while the total energy 
density and the total pressure are given by e = 2^ and p = 27^. 



APPENDIX D: RELATIVISTIC THERMODYNAMIC INTEGRALS AT ZERO 
TEMPERATURE, T=0 

At T = the mean occupation number is in good approximation described by the step 
function and the relativistic thermodynamic integrals Eqs. (A8) and (A9) becomes 
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The J n k integrals vanishes at T = because they are the product of temperature and 
the delta function or the X nk integrals. To calculate the integrals Eq. (Dl) we make the 
substitution p = msinhx. Then we have dp = m cosh x d;\\ If we set pf = msmhxf then 
the integrals Eq. (Dl) becomes 
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For massless particles Eq. (Dl) becomes 

7- 



-m n+2 / cosh n - 2fc X sinh 2(fc+1) X d X 
Jo 



(D3) 



4ttA 



(2k + l)\\J 



p n+1 dp 



(2A; + l)!!n + 2 (D4) 
The integrals Eq. (D3) are evaluated with the help of Eq. (2.413) of Ref. [25] and the 
properties of the hyperbolic functions. Here we list the results which are used in the text: 
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The Fermi momentum pj can be determined from the net charge density 
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from which Xf — sinh 1 (pf/m) 
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FIG. 1: The energy density and net baryon density dependence of X30 (in units of GeV 2 .fm~ 3 ) 
for the hadronic part of the equation of state. The energy density and net baryon densities are in 
units of the ground state densities, no = 0.16 fm~ 3 and Eq = 0.147 GeV fm~ 3 . 
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FIG. 2: The energy density and net baryon density dependence of J50 (in units of GeV 4 .fm~ 3 ) 
for the hadronic part of the equation of state. The energy density and net baryon densities are in 
units of the ground state densities, uq = 0.16 fm -3 and Eq = 0.147 GeV fm -3 . 
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